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ABSTRACT 

We present Fe Kev line profiles from and images of relativistic discs with finite thickness 
around a rotating black hole using a novel code. The line is thought to be produced 
by iron fluorescence of a relatively cold X-ray illuminated material in the innermost 
parts of the accretion disc and provides an excellent diagnostic of accretion flows 
in the vicinity of black holes. Previous studies have concentrated on the case of a 
thin, Keplerian accretion disc. This disc must become thicker and sub-Keplerian with 
increasing accretion rates. These can affect the line profiles and in turn can influence 
the estimation of the accretion disc and black hole parameters from the observed 
line profiles. We here embark on, for the first time, a fully relativistic computation 
which offers key insights into the effects of geometrical thickness and the sub-Keplerian 
orbital velocity on the line profiles. We include all relativistic effects such as frame- 
dragging, Doppler boost, time dilation, gravitational redshift and light bending. We 
find that the separation and the relative height between the blue and red peaks of the 
line profile diminish as the thickness of the disc increases. This code is also well-suited 
to produce accretion disc images. We calculate the redshift and flux images of the 
accretion disc and find that the observed image of the disc strongly depends on the 
inclination angle. The self-shadowing effect appears remarkable for a high inclination 
angle, and leads to the black hole shadow being completely hidden by the disc itself. 
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1 INTRODUCTION 

The fluorescent Ka iron emission line has been observed 
in several active galactic nuclei (AGN) with a broad and 
skewed line profile. The line is thought to be produced by 
iron fluorescence of a relatively cold X-ray illuminated ma- 
terial in the innermost parts of the accretion disc. Bear- 
ing in mind that the line is intrinsically narrow in the local 
rest-frame of the emitting material, and is transformed into 
broad, skewed profile by Doppler shifts and gravitational 
redshift effects, thus the line profile encodes the nature of 
the structure, geometry, and dynamics of the accretion flow 
in the immediate vicinity of the central black hole, as well as 
the geometry of the space-time, thereby providing key infor- 
mation on the location and kinematics of the cold material. 
Investigating these spectral features in X-ray luminous black 
hole systems opens a unique window allows us to probe the 
physics that occurs in the vicinity of a black hole, and pro- 
vides one way to test theory of strong field gravity. 
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Calculations of the line profiles emitted from an ac- 
cretion disc around a black hole have been performed by 
several authors. Theoretical Fe Ka line profiles from a 
thin disc around a S c hwarz s child b lack hole were calcu- 
lated by iFabian et all (|l989h . lLaorl (|l99ll l extended those 
to the extreme Kerr metric. These calculations are based 
on a geometrically t hin, optically thick accretion disc (here- 
after SSD, following Ishakura fe Sunvaevl l|l97St h on which 
the accreting material is considered to be in Keplerian or- 
bit around a central black hole. Further, the line emis- 
sivity is assumed to vary with r in power-law form. Ef- 
forts have been made later on to include various physi- 
cally plausible processes in the accretio n flow, such as spi- 
ral wave, disc warp, and disc thickness |Pariev fe Bromleyl 
19981 : lFukuell200dl : lHartnoll fe Blackmarj|2000l . l200ll . 120021 : 



Fukumura fe Tsurutdl2004h . as well as taking into consid- 



eration of the geo metry and the relative motion of the pri- 
mary X-ray source (Ruszkowski 2000; Dabrowski & Lascnby 
l200lj ; iLu fe Yull200ll : iNavakshin fe Kazanasll200ll ) towards 
a more realistic emissivity distribution. Some authors 
considered also the ionization effect, the emission from 
plunging region on the iron K line and reflection spec- 
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trum dMatt. Fabian fe Rossi Il993l; iRoss fc Fabianl 1 19931 : 



iMatt, Fabian fc Rosslll996l : iRevnolds fc Begermanlll997l ). 

In the calculations, two basic approaches have been 
used to map the disc into the sky plane. The first method 
follows photon trajectories starting from a given initial lo- 
cus of emission region in the local rest frame of the disc 
to t he observer at infinity. In this case a transfer func- 
tion llCunn ingham 1975; Laor 199l l: ISpeith. Riffert fc Ruder] 
1 19951 : 1 Wilms, speith fc Revnoldsf l998) is usually introduced 
as an integration kernel which includes all relativistic effects 
in line profile calculation. The integration for the line flux is 
then performed directly on the surface of the accretion disc. 
The t ransfer function was first introduced bv ICunninghaml 
l| 19751 ). who presented the numerical results for a grid of 
parameters aiming at estimating the relativistic effect on 
the continuum emission fro m SSD, and was re -fined and 
discussed in great detail by ISpeith et al l dl995h. The sec- 
ond method adopt a ray t r acing approach fD abrowsk i et al.l 
19971: iFanton et al.l Il997l: ICadez. Fanton fc Calvanil Il998l : 



Mtiller fc Camenzindll20o3 : ICadez fc Calvanill2005l ). Follow- 
ing the trajectories of photons from the sky plane to the 
accretion disc, in this method the image of the disc on 
the observer's sky is derived first and then the line flux is 
obtained by integrating over this image, weighted by the 
redshif t factor and the radial disc emissivity profile. Re- 
cently, [Beckwith~^^oi3 l|2004 ) developed a fast, accurate, 
high- resolution code which can b e used to generate high - 
resolution line profiles numerically. Beckwith fc Dong (|2005l ) 
extended it to include the contribution of higher-order pho- 
tons to the line profiles. But all of these approaches are 
restricted to SSD. 

On the other hand, direct imaging of accretion discs 
around a black hole is one of the most exciting areas of study 
to be considered in the future. The fact that a black hole 
posses an event horizon makes a black hole cast a shadow 
upon the background light with a size of roughly ten grav- 
itational radii that is due to the bending of light by the 
black hole, and this shadow is nearly independent of the spin 
or or ientation (Fal cke. Melia fc Agoll l2000; Zak harov et al.l 
l2005h . However, for a black hole embedded in an optically 
thick accretion flow the shape and position of t he shadow 
will be altered regardless of the black hole spin (|Takahashil 
12004 IWatarai et al.l [2005h . From an observational point of 
view, the highest angular resolution is obtained with Very 
Long Baseline Interferometry at millimetre wavelengths (so 
called mm-VLBI) with an angular resolution of a few ten 
micro-arcseconds. This corresponds to a spatial resolution 
of only a few ten gravitational radii for the nearby galaxies. 
Future Global mm-VLBI Array at short millimetre wave- 
lengths therefore should allow to map the direct vicinity of 
the Super Massive Black Holes (SMBH) such as Sgr A* and 
M87, and offers new p ossibilities to study the immediate 
environment of SMBH j Krichbau m et alj 12004: IShen et al.l 
120051 : iBroderick fc Loebll2006l : lYuan. Shen fc Huandl2006l l. 
In the X-ray band, the proposed Micro- Arcsecond X-ray In- 
terferometry Mission (MAXIM) aims to obtain submicroarc- 
second resolution X-ray i mages of n earby galactic nuclei 
l|Revnolds fc Nowakl I2OO3I : ICashl l2005h . At this resolution, 
one can capture the image of an event horizon of the cen- 
tral massive black hole in a nearby AGN. The combination 
of high resolution radio interferometry and interferometric 
X-ray spectroscopy would form an powerful tool to study 



SMBH and their environment with high accuracy and pro- 
vide unprecedented and direct constraints on the dynamics 
and geometry of the disc, as well as the geometry of the 
space-time. 

With the development of the observational techniques, 
high-quality observational data will eventually become avail- 
able. By fitting the data, one can in principle constrain the 
parameters of the accretion disc system, this will provide 
both a direct evidence for the existence of a black hole 
and a way of quantitative test general relativity in strong 
gravity. However, accurate quantitative analysis of obser- 
vational data requires a sophisticated model that treats 
all relativistic effects with a realistic accretion disc struc- 
ture. At present, such a complete model is still not avail- 
able. To our knowledge, SSD breaks down when the ac- 
cretion rate approaches the Eddington rate. At this limit 
the disc must become geom e trical l y thick and b e sub- 
Keplerain (lAbramowicz et alj Il988l : IWang fc Zlioul Il999l : 
IShadmehri fc Khaienabil l2005h. that is the so called slim 
disc. For a thick disc. iNandra et all fl995) pointed out it 
would be of lower density than a standard a— disc which 
would increase the ionization parameter (flux divided by 
density), thus leads to iron in the inner parts of the disc 
becoming fully ionized and no iron lines at all. For slim 
disc, this may not be the case, the broad, ionized Fe Kq 
line was discovered in some narrow-line Seyfert 1 galaxies 
(NLS1) dB allantync. Iwa sawa fc Fa bian 2001a; Boile r et al.l 
12001 120031 : iFabian et al l 120041 : iGallo et all l2007h . which 
have been thought to work with high accretion rates. So, slim 
disc has been received much more attention because of it can 
be basically used to account for spectral features in NLS1 
llMineshige et alfcood : IWang fc Neteerll2003l : IChen fc Wanel 
l2004 l. With the increasing evidence for ionized accretion 
discs in NLS1, the spectra and emission lines of slim discs 
need to be studied in more details. Motivated by the above 
considerations, a geometrically and optically thick accretion 
disc model is presented making an attempt at gaining an in- 
sight into the effects of disc geometry and dynamics on the 
lin e profiles and disc images. Following the idea presented 
bv lSpeith etHI (|l995h . we extend their method to the finite 
thick disc, and adopt elliptic integrals to improve the per- 
formance of the code which is much faster than the direct 
integral and widely used by many authors. 

The paper is organized as follows. In fj2] we summarize 
the assumptions behind our model, and present the basic 
equations relevant to our problem, while some more techni- 
cal aspects like the formulae of the integration for photon 
trajectories expressed in terms of the inverse Jacobian el- 
liptic functions are given in Appendix [A] We present our 
results in Sj3[ and summarize the conclusions and discussion 
in $H 



2 ASSUMPTIONS AND METHOD OF 
CALCULATION 

The aim of this paper is to consider how the accretion disc 
geometry and dynamics affect the iron Ka line profiles and 
disc images. To this end, the disc shape and structure must 
be determined first. To obtain a rigorous model, one should 
solve the disc structure equations numerically. However, this 
is beyond the scope of the current work. For simplicity, we 
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adopt a conical surface for the disc geometry. The thickness 
of the disc can be described by the half subtending angle 5 
(0 sC 5 < tt/4). When 6 = 0, the disc reduces to SSD. The 
complementary angle of 5 is denoted by i? c which is the an- 
gle between the symmetric axis of the system and the radial 
direction of the disc surface. The parameters of this model 
include: the radii of the emitting disc zone n n , r out , the spin 
of the black hole a, the inclination angle of the disc (i? ) and 
the disc surface angle (i? e ), the radial emissivity index p, the 
angular velocity index (see below) n, respectively. In addi- 
tion to all of those, the angular dependence of the emissivity 
also shall be given. 



2.1 Assumptions and basic equations 

The propagation of radiation from the disc around a Kerr 
black hole and the particle kinematics in the disc were stud- 
ied by many authors. We review properties of the Kerr met- 
ric and formulae for its particle orbits, and summarize here 
the basic equations relevant to this work. Throughout the 
paper we use units in which G = c = 1, where G is the 
Gravitational constant, c the speed of light. The background 
space-time geometry is described by Kerr metric. In Boyer- 
Lindquist coordinates, the Kerr metric is given by 

ds 2 = -e 2v dt 2 + e 2,p (d<f> - tudt) 2 + ^dr 2 + E(M 2 , (1) 



where 

2v 



ZA/A, e 2 ^ = sin" i?A/E, uj = 2Mar/A, 



E = r + a cos A = r + a 
A = (r 2 + a 2 ) 2 -a 2 A sin 2 0. 



2Mr, 



Here M, a are the black hole mass and specific angular mo- 
mentum, respectively. 

The general orbits of photons in the Kerr geome try can 
be ex pressed by a set of three constants of motion (|Carterl 
1 19681 ). Those are the energy at infinity E, the axial com- 
ponent of angular momentum EX, and carter's constant 
Q (=q 2 E 2 ). The 4-momentum of a geodesic has components 

Pm = (pt, Pt, pa, P4>) = (-E, ±EVR/A, ±EVQ, EX), (2) 
with 

R = r 4 + (a 2 — A 2 — q 2 ) r 2 + 2M [q 2 + (A - a) 2 ] r - a 2 q 2 
= q 2 + a 2 cos 2 § — A 2 cot 2 i? . 

From this, the equations of motion governing the orbital 
trajectory can be obtained. The technical details are given 
in Appendix lAl 

We assume that the disc is of a cone-shaped surface, 
axisymmetric, and lies on the equatorial plane of the black 
hole. Photons are emitted or reflected from the gas on the 
conical disc surface which moves along circular orbits. The 
radial drift of the gas on the disc surface is neglected. Thus, 
the 4- velocity field is chosen to be of the form 



u\d t , 0, 0, Vld^,) = (u, 0, 0, 



(3) 



where Q = vr/u* is the angular velocity of the emitting 
gas. The choice of Q must satisfy the causality condition. 
For sub-Keplerian v elocity, we adopt the modific ation of Q 
firstly introduced by iRuszkowski fc Fabianl (|200Ch 



n = 



7r/2 



l/r, 



i - 



tt/2 



l/r, 



(4) 



where ■& is the poloidal Boyer-Lindquist coordinate, fix — 
M 1/2 /(r 3/2 + aM l/2 ) is the Keplerian angular velocity and 
w is the angular velocity of the gravitational drag. It is easy 
to verify that Q. ^ fix- 

For describing physical processes near a Kerr black 
hole, Boyer-Lindquist coordinates, which are unphysi- 
cal in the ergosphere, are inconvenient. In order to 
make physics appear simple in their frames, the lo- 
cally nonrotating frames (LNRF) was introduced by 
iBardeen. Press fc Teukolskvl l|l972D . The relation between 
the local rest frame attached to the disc fluid and LNRF 
is given by a Lorentz transformation. In the LNRF, the az- 
imuthal component of 3-velocity reads 



,4sintf 

~e7a" 



(•>) 



The corresponding Lorentz factor 7 as measured by LNRF 
is defined as 7 = (1 — « 2 ) _1//2 . 

Due to relativistic effects, the photon frequency will 
shift from the emitted frequency v B to the observed one v a 
received by a rest observer with the hole at infinity. We in- 
troduce a g factor to describe the shift which is the ratio of 
observed frequency to emitted one: 



Vo/v c = (p ■ Uo)/(P ■ Uc 

e v (l-v 2 ) 1/2 /{l-VL\), 



(6) 



where p,u ,u c axe the 4-momentum of the photon, the 4- 
velocity of the observer and the emitter, respectively. 

The specific flux density F {v ) at frequency v Q as ob- 
served by an observer at infinity is defined as the sum of 
the observed specific intensities I (y ) from all parts of the 
accretion disc surface, 



F (v ) = / I (v )dQ. ohs , 



(7) 



where dfl ^ B is the element of the solid angle subtended by 
the image of the disc on the observer's sky. We do not con- 
sider the effect of higher order images of the disc in the 
following computations as their contribution is small due to 
most high order photons reintercept and be absorbed by the 
disc. 

Using the fact that / v 3 is invariant along the path of 
a photon, where v is the photon frequency measured by any 
local observer on the path , equation © can be rewritten as 



^0(^0) 



g I e (i/ e )dQ obs 



(8) 



I e (v e ) is the specific intensity measured by an observer coro- 
tating with the disc, and can be approximated by a <5- 
function, I' c (v' c ) = eS(v' c — v e ) where e is the emissivity per 
unit surface area. From well known transformation proper- 
ties of (^-functions we have 5(Vc — u Q ) = gS(v — gv c ), using 
this in equation (O, we obtain 



^0(^0)= / eg 4 5(v — gv c )dQ, oh 



(9) 



In order to calculate the integration over di} obs , we must 
first obtain the disc image or find the relation between the 
element of the solid angle and the disc linked by the null 
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geodesic. The apparent position of the disc image as seen by 
an observer is conveniently represented by two impact pa- 
rameters a and f3, measured relative to the direction to the 
centre of the black hole. The impact parameters a and j3 are, 
respectively, the displacements of the image in the directions 
of perpendicular and parallel to the projection of the black 
hole s pin. They are related to two constants of motion A and 
q by (|Cunningham fc Bardeenlll973l ; ICunninghamlll975h 

a = —A/ sin # , (3 = ± (q 2 + a 2 cos 2 i9 — A 2 cot 2 i9 ) , (10) 

where i9 is the angle between the observer and the rotation 
axis of the black hole (i.e. the inclination angle). The element 
of solid angle seen by the observer is then 

dad/3 1 d(a,(3) d(X,q) 



dQ a 



rl r% d(X,q) d(r,g) 



dr dg 



8(X,q) 



To/3 sin # d(r,g) 



dr dg, 



(11) 



where r is the distance from the observer to the black hole. 

Substituting equation into equation (|9]) gives the 
desired result: 



F {u ) = 



r 1/3 sin tf D 



eg 4 S{v - gv c ) d }, X,q l dr dg. (12) 
d(r,g) 



To perform the integration, the form of the disc emis- 
sivity in the integrand also needs to be given. In general, it 
can be a function of the radius, r a and polar angle, n c , of 
an emitted photon with the surface normal of the disc in 
the rest frame of the emitting gas. This angle is determined 
by taking the dot-products of the photon four-momentum p 
with the surface normal n. The surface normal in the rest 
frame is 



n = Y,- 1/2 d/d$. 

By definition, we get 

p ■ n p ■ n p ■ u 

P ■ U e p ■ Uo p ■ U e 



(13) 



cos(n c ) = 



= gVO/Vz 



e"(l - u 2 ) 1/2 e 1/2 E- 1/2 /(l - SIX). 



(14) 



If the emission is isotropic in the rest frame, we do not need 
to know n e . More generally, we take the form 



e(r c , Me) = e(r e )/(/x e ), 



(15) 



where /i c is the cosine of the emitted angle (cosn e ). And 
the radial emissivity is assumed to vary as a power law with 
emissivity index p: 

e(r c ) oc r~ p . (16) 

We consider three p ossible cases for th e angu lar dependence 
of the emissivity (Bec kwith fc Done] [2005): (1) isotropic 
emission, /(/ie) = 1; (2) limb darkening law /(/i c ) oc (1 + 
2.06/i c ); (3) limb brightening law /(/i c ) oc l/fi c . 



by n the radius of each such emitting ring. For each ring 
there is a family of null geodesic along which radiation flows 
to a distance observer at polar angle $ from the disc's axis. 
As far as the iron Ka emission line is concerned, for a given 
observed frequency v a the null geodesic in this family can be 
picked out if it exists. So, the weighted contribution of this 
ring to the line flux can be determined. The total observed 
flux can be obtained by summing over all emitting rings. 
Changing the observed frequency, then the line profiles will 
be obtained. 

This family of null geodesic for each ring can be used to 
map the accretion disc onto sky plane, that is, disc imaging. 
A geodesic in this family connects an emitting region in the 
ring to the distant observer. The constants of motion A and 
q of this geodesic can be used to determine the apparent 
position of the emitting point on the sky plane using the 
corresponding two impact parameters a and 13. Different 
geodesic is associated with different point. Using geometric 
optics, one determines the appearance of the ring from this 
family of geodesic, then in this way images (at infinity) of 
the accretion disc are obtained. 

The main numerical procedures for computing the line 
profiles are as follows: 

(i) Specify the relevant disc system parameters: 
fin, r out , a,p, n, i9 , i9 and the angular emissivity. 

(ii) The disc surface is modeled as a series of rings with 
the radii n and weights uii which calculated u sing an algo- 
rithm due to Rybicki G. B. (|Press et al.lll992j '). 

(iii) For a given couple (n, g) of a ring, the two constants 
of motion A and q are determined if they exist. This is done 
in the following way: the value of A is obtained by the an- 
other form of equation © 



a = Ui-*L =1 

Q V jg) n 



1 - 



e »(l_^)i/a 



(17) 



the value of q is determined for solving photon trajectory 
equation (|A1[) . Then the contribution of this ring on the 
flux for given frequency v with respect to g is estimated. 

(iv) For a given g, the integration over r of the equa- 
tion (|12[) can be replaced by a sum over all the emitting 
rings 



Fo{Vo) 



4 

qev Q 



8(X,q) 



^ rlvt (3 sintf d{r,g) 



(18) 



The Jacobian [<9(A, q) /d(r, g)] in the above formula was eval- 
uated by the finite difference scheme. From above formula, 
one determines the line flux at frequency v a from the disc, 
(v) Varying g, the above steps are repeated. 

The observed line profile as a function of frequency v 
is finally obtained in this way. 



2.2 Method of calculation 

With all of the preparation described in the previous section, 
we now turn to how to calculate the line profiles and the 
disc images numerically. We divide the disc into a number 
of arbitrarily narrow rings, and emission from each ring is 
calculated by considering its axisymmetry. We shall denote 



3 RESULTS 

The model and the computer code described in this paper 
suitable for disc inner edge located at any n n ^ r ms , where 
r ms is the radius of the marginally stable orbit. For sim- 
plicity, in all plots presented in the paper, we assume that 
fin = f ms . We have taken a disc from r ms to r max = 20r g (fo- 
cus on strong gravitational effects) for Kerr metric case 
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Figure 1. The relativistic line profiles computed by our code 
both for the Schwarzschild (a = 0) and maximal Kerr (a = 0.998) 
cases for 1? = 30° (left) and 85° (right). The disc zone is from 
r ms to r ou t = 20r g and located at equatorial plane, where r g is 
the gravitational radius. Upper panel: The Schwarzschild metric 
for e(r c ) oc r^ 3 and /(/J c ) = 1. Lower panel: The maximal Kerr 
metric for e(r c ) oc r~r 3 and f(fie) oc (1 + 2.06/i c ). The flux in all 
cases is given using the same arbitrary units, and all our results 
are unsmoothed. 



Figure 3. The relativistic line profiles as a function of the angular 
velocity represented by the parameter n in equation 10}: Q = 



w for n 



1,2,3,4 (from top 



to bottom at the redshift peak) for a maximal Kerr black hole 
with the disc extending from 1.235 — 20r g and i9 = 30°,i9 c = 
60° (upper panel), 70° (lower panel). The emissivity law is the 
same as in Fig. [2] 
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Figure 2. The relativistic line profiles as a function of the 
disc thickness for (from bottom to top at the red peak) -d c = 
50° , 60° , 70° , 80° and 90° for a maximal Kerr black hole with the 
disc extending from 1.235 — 20r g . The observer inclination equals 

30° and angular velocity takes the form: Q = (^75) ^ K ~^ 

cj, here n is set to 2 (upper panel) and 3 (lower 

panel). The emissivity law is taken the forms e(r c ) oc r^T 3 and 
/(/ie) oc (1 + 2.06/ie). 



treated spin of the black hole as a free parameter for dif- 
ferent observed inclinations and disc thickness. Due to its 
astrophysical importance, we choose the iron fluorescence 
line at 6.4 keV in what follows. 



3.1 Relativistic emission line profiles 

The numerical code discussed in previous section was used 
to model emission line profiles for different model parame- 
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Figure 4. The relativistic line profiles as a function of the ob- 
served inclinations for i9 = 15° , 30° , 45° , 60° and 75° for a max- 
imal Kerr black hole with the disc extending from 1.235 — 20r g 
and j? e = 60° (upper panel), i9 e = 70° (lower panel). The index 
n is set to 3 and the emissivity law is the same as in Fig. |2] 



ters. To test the performance of our code, we first compared 
the line profiles generated by our code when the disc is re- 
duced to SSD to those g enerated by the code described in 
iBeckwith fc Do ne (2004), and found that the overall match 
is fairly good, especially for the Schwarzschild metric case. 
FigQJsho ws the results for paramet ers identical to those pre- 
sented by IBeckwith fc Pond (|2004h . 

The dependence of the line profiles on the disc thickness 
is shown in Fig[2] The angular velocity and the emissivity 

law take the forms: Q — ^'^75^ 

e(r c ) oc r^ 3 and f(p c ) oc (1 + 2.06/i c ). Two cases for 
n — 2 (top panel) and n = 3 (bottom panel) are presented 
in this figure. It is explicit that the separation and relative 
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Figure 5. The relativistic line profiles generated by our model 
with (a) e(r c ) oc r<r 3 , f (/i c ) oc (1 + 2.06/i e ) (black line), (b) 
e (r e ) oc r~ 3 , / (jtt e ) = 1 (red line), (c) e (r e ) oc r~ s , f (fj, e ) oc /x" 1 
(green line), for a maximal Kerr black hole with the disc extending 
from 1.235 - 20r g and i? Q = 30°,i3 c = 70°. The sub-Keplerian 
angular velocity is the same as in Fig. [2] All profiles are scaled to 
unity for better comparison in this case. 
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Figure 6. Comparison of the relativistic line profiles generated 
by our model with different spins a = 0, 0.5, 0.998. The emission 
line region is from r ms to 20r g and the angular velocity and the 
emissivity law are the same as in Fig. [4] The angles are marked 
in each figure. 

height between the blue and red peaks diminish as disc thick- 
ness increases because the disc becomes more sub-Keplerian. 
This effect is also clearly illustrated in FigO The index n in 
Q describes the deviation of the angular velocity from Keple- 
rian one. As the deviation from Keplerian velocity increases 
the height of the blue peak of the line decreases significantly. 

FigSl compares the line profiles at different viewing 
angles # = 15° , 30° , 45° , 60° and 75° for a maximal Kerr 
black hole with the disc extending from r ms to 20r g and 
i?o = 60° , 70° . At high inclinations the self-shadowing effect 
has been taken into account. Due to gravitational lensing 
(light bending) effect, there is still a substantial fraction of 
light that can reach the observer at infinity. 

We also calculated the effects of emissivity on the form 
of the relativistic line profile. The radial emissivity is taken 
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Figure 7. As in Fig. [6] but with the angle i9 c = 70°. 
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Figure 8. As in Fig.|6]but with the angles i9 Q = 85° and tf c = 70° 
(upper panel), 60° (lower panel). 

a power law with the index p, which determines the relative 
contribution from different radii of the disc. Here we focus 
on the influence of anisotropic emission on the line profile. 
Different angular emissivity laws have striking effects on the 
line profile, which we illustrate in Fig[5]for a maximal Kerr 
geometry with the disc extending from r ms to 20r g and i9 = 
30° , i9 e = 70° . The angular emissivity takes one of the three 
forms: (a) e (r c ) oc r<T 3 , / (jj, e ) oc (l + 2.06/i c ), (b) e (r c ) oc 
nT 3 , /(^c) = 1, (c) e(r c ) oc r~ 3 , / (^ c ) oc fj,~ l . From the 
figure one can see the relative height of the blue wing changes 
a lot for different angular emissivity laws, anti-correlated 
with the slope of the red wing. 

The line profiles as a function of the black hole spin 
are also demonstrated. For a low or intermediate inclination 
angle the line profiles are shown in Figs [S] and [7] Note that, 
the red wings change significantly whereas the blue peaks 
almost are not affected by the spin. At high inclinations, 
the effect of the self-shadowing dramatically alter the line 
profile for a thick disc. The results are illustrated in Figs 
[8] and [9] with angular emissivity f(fj, e ) oc (1 + 2.06^i e ) and 
f(^e) = 1, respectively. For § c — 60°, the line profiles are 
almost the same, this implies that the line emission from 
the inner parts of the disc are completely obscured by the 



Line profiles from thick discs 7 




0.10 
0.05 
0.00 



(i3 =S5°, i9 e = 60°) 

Q = 

a = 0.5 

a = 0.99S 



0.6 0.3 

g-factor 




Figure 9. As in Fig. [8] but with no angular dependence of emis- 
sivity e (r e ) oc r^ 3 , / (fj, e ) = 1. 



outer parts of the disc. At high viewing angles, the impact 
of angular emissivity law on the relativistic line profiles are 
also striking. 



3.2 Accretion disc images 

We present in Figs [TU] and [TT]the redshift and flux images of 
the accretion disc and black hole shadows on the (a, (5) plane 
for an extreme Kerr black hole, for $ = 5°, 30°, 55°, 80°, 
and i? c = 60° , 70° . Redshift images are colored by the as- 
sociated values of g as measured by the infinity observer, 
which is defined by the scale at the top of each image. Flux 
images are colored by 10 4 eg 4 , again with the scale defined at 
the top of each image. The parameters n and p are both set 
to 3 and f(fi e ) oc (1 + 2.06/i c ). The images are distorted by 
the combined action of Doppler effects, gravitational redshift 
and light bending in the vicinity of the black hole. Note that 
at small inclination angle ($ = 5°), the observed radiation 
is all redshifted, and therefore the emission line profiles will 
have a net redshift. On the other hand, at an intermediate 
inclination angle ($ = 30°, 55°), the innermost part of the 
disc is notably redshifted, whereas the observed radiation 
from the approaching side is remarkably enhanced by the 
Doppler boost. Moreover, the light ray emitted by the far 
side of the disc is bent by the gravity of the black hole, re- 
sulting in the vertical asymmetry of the image, as if it were 
bent toward the observer. Note also that the self-shadowing 
effect is remarkable at a high inclination angle (i9 = 80°), 
and therefore the black hole shadow in this case does not 
appear at all. The shape, size and position of the black hole 
shadows are also affected b y the self-shado wing, which is 
different from those of SSD (|Takahashill2004) . 



any value of the black hole spin, the different viewing an- 
gle, the disc inner radius (ri n r ms ), and the disc thick- 
ness (<5 ^ tt/4). It also allows the user to choose one of the 
three types of angular emissivity laws: isotropic emission, 
limb-darkening or limb brightening laws. 

We show that the separation and the relative height 
between the blue and red peaks of the line profiles diminish 
as the thickness of the disc increases because of the sub- 
Keplerian motion. The angular emissivity form has also a 
significant influence on the line profile. The results of one 
peak line profile present in intermediate viewing angle in 
our model is different from those in SSD for low viewing 
angle. To see the self-shadowing effect more clearly, images 
of the disc and the black hole shadows are also presented 
in this paper. The self-shadowing effect is very important 
for high inclination angle. Future X-ray observations of high 
state accreting systems such as narrow line Seyfert 1 galaxies 
will be important to test whether the disc in these systems 
are indeed thick. 

Here we just present a simple disc model with a con- 
ical surface aimed at getting an insight into the effects of 
geometric and dynamic influence on the line profiles and 
disc images. For a non-equatorial disc, we consider the self- 
shadowing and sub-Keplerian effects on them, as well as 
the contribution of the light with carter constant Q < 
which different from those for equatorial disc. For simplic- 
ity, in this paper we neglected the influence of radial drift 
of a flow on the line profile. Other limitations include the 
thickness of the disc may vary with radius, it probably also 
has a substantial warp, and the effects of photoionization 
of the surface layers of the accretion disc on the emis- 
sion lines are not taken into account. X-ray reflection by 
photoio nized accretion discs has been investigated in some 
detail dRoss fc Fabianl Il993l. 120051; iRoss. Fabian fc Youngl 
ll999l ; lBallantvne. Ross fc Fabianll2001br ). The ionization pa- 
rameter has clearly a large effect on emissions lines. Ev- 
idence for reflection by ionized accretion discs in NLS1 
has been accumulated in the l i terature in recent years 
(see e.g. iBallantvne et al.ll2001al: iBoller et all |2002| . 120031 ; 
iFabian et alj 12004 iGallo et al.ll2007l ). Furthermore, the ra- 
dial drift of the flow for a sub-Keplerian disc may has sig- 
nificantly influence on the line profiles. A more realistic disc 
model should take into account both the sub-Keplerian and 
radial velocity effects on the line profiles. This effect will be 
investigated in the near future. 
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4 SUMMARY 

We have developed a computer code both to calculate the 
line profiles of a relativistic thick accretion disc around a 
black hole and to generate the images of accretion discs. 
The code includes all relativistic effects. It also includes the 
effect of self-shadowing of the disc, i.e. the outer disc blocks 
the emission from the inner region. The code can handle 



8 Wu and Wang 




•» -m p- »p w -bp -ip p 



It M -» -19 9 W N -» -IP P 19 » 



in -hi su uli iw as mi hi hi m «i 




-J9 -« 9 IP 90 



-» -19 9 ■# *9 -» -«P P 19 ,M 



Figure 10. Rcdshift images (top row) and flux images (bottom row) of the accretion disc on the (a, (3) plane for a extreme Kerr black 
hole. From left- to right-hand side: i9 = 5°, 30°, 55°, 80° and i9 e = 60°. Redshift images are colored by the associated values of g as 
measured by the infinity observer. Flux images are colored by 10 4 eg 4 . The parameters n and p are both set to 3 and f(n c ) oc (l + 2.06/x c ). 
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APPENDIX A: INTEGRATION OF PHOTON ORBITS IN A KERR SPACETIME 

The general orbit of a photon (indeed for any particle) in a Kerr spacetime is described by three constants of motion l)Carterl 
1968): the energy-at-infinity E, the angular momentum about the axis of the black hole L z , and carter's constant Q . Let 
us write L z — XE and Q — QE 2 . Then, the equations governing the orbital trajectories are separable. Since the system 
is stationary and axisymmetric, only the motions in the r- and ^-dire ctions are required in the calculation of th e radiation 
spectrum from the disc. The motion in the r — 1} plane is governed by (jBardeen et al.lll972l ; IChandrasekharll 19831 ) 



dr 



where 
R(r) 

0(0) 



dO 



4 , / 2 

r + [a 



' + a 2 cos 2 ' 



Q) r 2 + 2[Q+(\-af]r-a 2 Q 



A cot 1? , 



and r e and # e are the starting values of r and 0. 
Define fi = cos then equation (| Al|l becomes 



dr 



= ± 



dfi 



where fi e — cos$ e and 

@fj,(p) = Q + (a 2 — A 2 — Q)fJ? — a 2 fi 4 — a 2 (/j 2 + fi 2 
and n± are defined by 
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(A6) 



2 2 2 2 2/2 

For Q > 0, both fi + and fi_ are non-negative, when Q < 0, p_ is less then zero. Note, fJ.+fJ,- — Q/a . 

For a photon emitted by the disc, the integral over fi can be worked out with the inverse Jacobian elliptic integral 
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(A7) 



where ^ fi < fi+ for fi 2 _ > and \J —\x % _ K fx < fl+ for /i"! < 0. 

The integral over r can also be worked out with inverse Jacobian elliptic integrals. To do so, we need to find out the four 
roots of R(r) — 0. R(r) = may has four real roots; two real roots and two complex roots or four complex roots. We consider 
the three cases separately: 

Case A. R(r) = has four real roots Let us denote the four roots by ri, r 2 , 7-3 and 7-4 in the descending order. 
The integral over r can be worked out by the following integration 



dr 

Vr 



dr 



n y/(r - ri)(r - r 2 )(r - r s )(r - r 4 ) 



where 



7714 = 



\J (ri - r 3 )(r 2 - rt) 



r4.)(r 2 - 



rz) 



(n — r3)(r2 — r&) 



(r 2 — ri)(r — ri) 



(n — n)(r — r 2 ) 



Til 4 



(A8) 



(A9) 



When n =V2, the integral over r can be expressed in terms of a logarithm ,which is of no practical interest for this paper. 

Case B. R(r) — has two complex roots and two real roots Let us assume that ri and r 2 are complex, 7-3 and 7-4 are 
real and r 3 > r& . The physically allowed region for photons is given by r > 7-3 . If we write n and r 2 in the form n =u + iv 
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and n = u — iv , the integral over r can be worked out with the following integration 
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y/[{r — u) 2 + v 2 ](r — rs)(r — r*j 
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(All) 



K(vti2) is the com plete elliptic integral of the first kind. When r = oo, equation (IA10|) is equivalent to equation (28) of 
Cadcz et al.l (jl998h . It can be directly verified that 



(A12) 

Case C: R = G has four complex roots Let us denote the four roots by n, r2, r 3 and r 4 . n = rj, 7*3 =r%, where * stands 
for complex conjugate. If we set B = |ri + r-2 1 , C = riT2, -D = r3r 4 , which leads to i?(r) = (r 2 — Br + C)(r 2 + Br + D). In 
this case, in order to express the quartic in terms of r , we make the substitution 

r — A 2 
r - Ai 

where Ai, A 2 (Ai < A2) is the two real roots of equation 
2BA 2 + 2(D - C)A - B(D + C) = 0. 
This yields 

(pis 2 + gi)(p 2 s 2 + g 2 ) _ pip 2 (s 2 +p 2 )(s 2 + g 2 ) 



R(r)~- 
where 



(8-1)* 



Pi = A 2 — BAi + C, 
p 2 = A? + BAi + D, 



(s-iy 

qi = Af - BA 2 + C*, 
g 2 = Al + BA 2 + D, 



(A13) 



p = max{g 1 /pi,g 2 /p 2 }, g = min{gi/pi, g 2 /p 2 }. 



From equation (|A13|) the integral over r can be calculated by means of the substitution s = through the following 

formula 
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where 



